
% Annual
m=nmodel.data;
d=dat.data;
di=d-m;
load cmap_insitu
close all
gccolor('mp',dat.londata,dat.latdata,m','on',0,1000,['EKE - Model - Annual - ',year,' (cm-2.s-2)'])
colormap(cmap_insitu)
printgc('-f1','-depsc2',['./PLOT_VALIDATION/EKE/EKE_model_',parameter,'_',bbio,'_y',year,'_ANNUAL.eps'])
close all
gccolor('mp',dat.londata,dat.latdata,d','on',0,1000,['EKE - Data - Annual - ',year,' (cm-2.s-2)'])
colormap(cmap_insitu)
printgc('-f1','-depsc2',['./PLOT_VALIDATION/EKE/EKE_data_',parameter,'_',bbio,'_y',year,'_ANNUAL.eps'])
load cmap_whitemiddle
close all
gccolor('mp',dat.londata,dat.latdata,di','on',-1000,1000,['EKE - (data-model) - Annual - ',year,' (cm-2.s-2)'])
colormap(cmap_whitemiddle)
printgc('-f1','-depsc2',['./PLOT_VALIDATION/EKE/EKE_diff_',parameter,'_',bbio,'_y',year,'_ANNUAL.eps'])
close all

